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Abstract 

The B — > wlv weak decay form factors are calculated via light-cone sum rules within the frame- 
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I. INTRODUCTION 



It is well-known that the weak B — ► tt transition is important for the extraction of 
the Cabbibo-Kobayashi-Maskawa matrix element V u b from semileptonic decays, and for the 
measurement of CP violation from non-leptonic decays. In order to achieve these goals 
advances have to be made both on the experimental and theoretical sides. On the theoretical 
side an important task is to reduce the uncertainties in the calculations of the relevant 
hadronic matrix element, represented by 



(w(p)\urb\B(P)) = f + (q 2 



(P + pY 



9 H 

q 



m B ~ m l u 



where q = P — p is the momentum transfer to the leptons. Up to now, the form factors 
have mostly been calculated via light-cone (LC) sum rules 0] in full QCD. Instead, here we 
will work in the context of the heavy quark effective theory (HQET) || . The reasons are as 
follows: 

(i) Unlike for heavy-to-heavy transitions an application of HQET methods does not sig- 
nificantly simplify the analysis for heavy-to-light transitions. Nevertheless, the systematic 
nature of HQET allows one to identify and estimate uncertainties in the heavy-to-light 
transitions more easily. 

(ii) For an analysis involving the B meson HQET is the correct approximation method 
of QCD. In certain parts of the phase space of the semileptonic decay, the pion is not very 
energetic so that HQET may still be valid. 

(iii) The B — > tttt amplitude has been shown to be factorizable at the leading order of 
the l/m b expansion [Q]. Knowledge of the form factors based on HQET is needed so as to 
consistently apply the factorization approach. 

(iv) Furthermore, it is still controversial whether the time-like transition B — > tt at large 
recoil is governed by perturbative or non-perturbative QCD Q. If non-perturbative effects 
dominate, then the use of HQET for the B — > tt transitions can be fully justified. 

(v) In addition, the results of HQET can also be applied to D — ► tx{K) transitions. In 
HQET, a model-independent analysis of B(D) — > tt transitions to order 1/mg has been 
done in Ref. g i and for B(D) -> p in Ref. [0. 

HQET simplifies the analysis by introducing a set of universal functions. However, in 
order to obtain information on the universal functions themselves, some nonperturbative 
techniques, such as light-cone sum rules or lattice simulations, must be used in addition. 
Here we adopt the LC sum rule method || [| which is suitable for the calculation of form 
factors when light energetic hadrons are involved. To our knowledge, the LC-HQET sum 
rules were first applied in Ref. |l(J. In Ref. [llj|, they were used at the leading order of 
HQET for B — > tt transitions. However, the results in this paper differ from those in p~l] . 



Note that we distinguish between LC-HQET sum rules and LC-QCD sum rules with the 



1 / mQ expansion [12| . The main reason is that the ways to include radiative corrections are 
different. Another reason is that for H — > tt transitions, there is a subtle difference between 
the two types of sum rules, as will be discussed in the paper. 

In this paper we apply LC-HQET sum rules to the H — > tt transition to order I/tuq. In 
the next section we give a brief review on the application of HQET to the decay B — > ttIv. In 
section [CTTJ, the leading and next-to-leading order universal functions are calculated by using 
LC-HQET sum rules. In Section |V| we present our numerical results. Section [V] contains a 
summary and discussion, and a comparison with other approaches. 
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II. THE HEAVY QUARK EXPANSION 



In HQET, the velocity of heavy quark Q, v, is a well defined quantity The heavy quark 
field can be represented by the velocity-dependent field, 



h v (x) = exp(irriQ v-x) P + Q(x) , 



(2) 



where P+ 



projects onto the upper component of the heavy quark field Q(x). To 



the order 1/mg, the effective Lagrangian is given by [13 

1 



£hqet = h v iv-Dh v + 



2mr 



Okin + CmagO) O 



mag 



(3) 



where the gauge- covariant derivative = <9 M — ig s T a A a generates the residual momentum 



and 



Okm = K (iD) 2 h v 



D - — h n 

^mag ll 'v u ap y ~ r ,L v • 



(4) 



Okin describes the kinetic energy of the heavy quark in the hadron, and O mag the heavy quark 
chromomagnetic energy. The equation of motion, i v ■ D h v — 0, is exactly satisfied. The 
higher dimension operators are treated as perturbative power corrections. In the leading 
logarithmic approximation, the renormalization factor C mag (/x) is given by 



CmagO) = r(/i)" 



r(//) 



(5) 



a(m Q ) 

where rif is the number of quarks lighter than the heavy quark Q. 

In the framework of HQET, it is convenient to work in the matrix representation for the 
description of the hadrons, in which wave functions of heavy hadron are only dependent on 
the heavy quark symmetry and their Lorentz transformation properties. The ground-state 
pseudo-scalar and vector heavy mesons are described by the so-called spin wave function 



M(v) = P 4 



"75 



for J p 
for J p 



0" 

r 



(6) 



with e M being the polarization vector. The form factors are considered as functions of the 
kinematic variable 



v ■ p 



2m 



H 



(7) 



Using the mass-independent normalization of the heavy meson state \H(v)) = m£- ' \H(P)), 
the form factors can be re-defined as 



(7r(p)| «7m Q \H(v)) = 2 if\(v ■ p) v, +f 2 (v p) %\ , (8) 
where the dimensionless variable is p 1 " 



The relation between the form factors in 



Eqs. ([J) and (||) is given by 



v ■ p 



fo 



2 



f2(v-p) fl(v-p) 



v ■ p 



m 



H 



ml 



fi(v-p) + h{v-p) 



v ■ p 



(9) 

fi(v-p)+p 2 f 2 (vp)}}.(10) 
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The form factors can be systematically analyzed order by order in terms of powers of 
1/ttiq. Let us begin with the expansion of the heavy-light vector current, 



(rr(p)\qYQ\H(v)) = ^Q^) <7r(p)| Ji\H{v)) 



+ E 



2m, 



-B^){v(p)\0 3 \H(v)) +0(- 



mf 



(11) 



In the limit of massless light quarks, a convenient basis of the above operators is WA 

qv^h v , 

q (—iv-D) ^h v , 
q (—iv-D) v^hy , 
q(-iD»)h v . 



Jl 


= qi^hy , 




o 1 


= q^iphy, 




o 2 


= qv^ ip h v , 


o 5 


3 


= qiD^hy , 


o 6 



(12) 



The corresponding Wilson coefficients are given by |T^, |I~5 

1 



B 2 (p) = - B 3 (y) = C 2 (p) = , 



34 



16 



■ r 



r 



28 2 uu -1 

— r H r 

27 27 



r lnr 



S 6 (/i) 



(13) 



10 

27' 27' Y + T 
20 

~9~' ~ UKr ~' " 3 

At the leading order of the 1 / rriQ expansion, the matrix element of the relevant current 
qT h v can be written as 



-2r 2 
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(tt(p) I g r ^ = -Tr{ 75 [ L> • p, p) + #L b (i> • p, p) ] V M(v) } 



(14) 



where the universal functions L Q (v ■ p, p) (a = a, b) depend on the kinematic variable v ■ p, 
but not on the heavy quark mass tuq. 

At the next-to- leading order in the heavy-quark expansion, the 1 / rng corrections coming 
from both the effective current and the effective Lagrangian of HQET will appear as follows. 
Matrix elements of the operators 0\, 2 , and O3 in the effective current can be expressed 
by the generic structure 

(7r(p)|g(riD) a 

= -Tr{[{F lVtl + F 2 p fl + F 3lfl ) l5 + {FtVp + F s Pr + F 6 iJi s p]rM(v)}, (15) 

where the universal functions Fi(v-p,p) (i = 1, • • • , 6) are also mQ-independent. Matrix 
elements of the operators 4 , 5 , and 6 are not independent, and can be obtained from the 
above structures Eqs. (|T4T) and fll5|) . Additionally, corrections coming from insertions of the 
operators Oki n and O mag into matrix elements of the leading-order currents can be described 
by six additional universal functions K a (v ■ p, p) (a = a, b) and Si(v ■ p, p) (i = 1, 
which are defined by matrix elements of the time-ordered products 



7T 



(p)\ij d A y T{qT h v (0), O kin (y)} \H(v)) = -Tr{ 7s ( K a + pK b ) T M(v) } 
(p)\ i f d 4 yT{qVhy{Q), O mag (y)} \H(v)) 



■A), 

(16) 



-Tr 



{iSlPa^/3 + S 2 CT af3 ) 75 + (iS 3 PoTf/3 + S 4 (Tap) 75 $ 



TP + a af3 M(v 



(17) 
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The form factors in Eq. (|8]) can then be expanded as follows after doing the appropriate 
traces IHI, 



fi — Ci L a 



f 2 = CiL b + 



1 



2m ( _ 



2mr 



G x ^ + {K a + C mag S a ) -B 4 {A-vp)L a - B 5 {A-vp) 



x(L a + L b ) + B % [Tl-K{L a + L b )\ | + 0( 
d \F l h +{K b + C^ g S b )] -B 4 (A-vp)L b 



2 > ' 



Tl + v ■ p {L a + L b ) + 0{ 



m 



2 / ' 



(18) 



(19) 



where [i = 1, 3) and S a are defined by 

^ = F 1 + 2F 3 -p 2 F 5 , 

T\ = F 2 + F 4 + 2F 5 - 4F 6 
S a = -2S 1 + 6S 2 + 2p 2 S 3 



F b = F 2 + F 5 — F e , 
S b = 2S 1 - 2S 3 + QS A 



(20) 



It should be mentioned that consequences of the heavy quark symmetry and the equations 
of motion for the heavy and light quark fields make only two (F 5 and _F 6 ) of the six form 
factors Fj's independent 0. This means that one can re-express T l a by L a , L b and F 6 as 
follows 



J a 
r b 



-(A - 2v ■ p)L a + v ■ pp 2 L b + 4F 6 
-vpL a -AL b -AF 6 , 



^ = v P L a + AL b + 2F 6 , 
Tl = -vpL a -AL b -2F 6 . 



(21) 



Note that between the two independent universal functions from the effective current cor- 
rections only Fq(v-p,^l) is relevant to the H — > 7T matrix element in Eq. ([I]) 



(22) 



III. THE LIGHT-CONE HQET SUM RULES 

Our aim is to calculate the independent form factors given in the last section by LC-HQET 
sum rules. In subsection A, the leading universal functions L a (a = a, b) are calculated. 
The relevant sub-leading universal functions, F 6 and K a + C mag S a , will be calculated in 
subsections B and C, respectively. 

A. Leading order 

Let us begin with the following 2-point vacuum-pion correlation function as depicted in 
Fig. 0: 



F^X^u) = i d 4 xe 



Ak-x 



7T 



T{u(x) lfM h v (x), h v (0)i l5 d(0)}\0) 



F a( X ^)v^ + F b (X,u;)p^ 



(23) 
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FIG. 1: The diagrammatic representation of the 2-point vacuum-pion correlation function Eq. n23), 
where the heavy solid line represents the heavy quark. 

where A = 2v -p and u — 2v • k. On the one hand, by inserting a complete set of intermediate 
states with the same quantum numbers as the H meson between the current in the above 
vacuum-pion correlation function, we obtain the hadronic representation 

F* adr (A,£) = 2F L ^ A l + resonances, (24) 

2A — LU 

where u — 2v ■ k with k being k + p. A is the heavy meson mass as defined in the HQET. 
Note that it is u) that is the relevant energy of the quark system. The decay constant F is 
defined as 

(H(v)\h v Td\0) = ^F(ji)Tr[TM(v)] , (25) 

with T being an arbitrary Dirac matrix. 

On the other hand the correlation function can be calculated by expanding the T-product 
of the currents near the light-cone at x 2 = 0. In our calculations, the chiral limit p 2 = m 2 = 
is taken. By adopting the fixed-point gauge in which 

x^{x) = , (26) 

the following complete heavy quark propagator equals the free one, 

r°° 1 + is 

(0\T{h v (x) h v (0)}\0) = / dt8(x-vt)-^. (27) 

J £ 

The matrix elements of the nonlocal vacuum-to-pion operators can be parameterized by pion 
wave functions on the light-cone. The LC wave functions are classified according to their 
twist and the number of partons they describe. Here, let WJ-'in) denote the n-particle light- 
cone wave function. One can easily find that at the leading order of HQET only the W^F^s 
are relevant. Up to twist four, the matrix elements are parameterized as follows 

(n(p)\u(x) W d(0)\0) = -%U Vil [\ue iupx {cp 7r (u) + x 2 g 1 (u)) 

2 

X Pp\ I .1,. ^u V x l 



-fnix^-^^-) j due 1 - //.,(//) 



px ' JO 
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(7r(p)|«(a;)i75d(0)|0) = f du e lu *») , 

J 

(n(p)\u(x)a^ l5 d(0)\0) = i ^ (p,x u - p u x,) £ du e iu ^<p a (u) , (28) 

with p n = m 2 / ' {m u +m d ). Here, (p n is the leading twist- 2 wave function, g 1 and g 2 are twist-4 
wave functions and tp p , ip a twist-3 wave functions. In the fixed-point gauge, the operator 
n G = P exp {ig s Jq da x II A fl (ax)^ is unity. Considering the identity 

l^lu = -i<V + SV > ( 29 ) 
one can obtain F^ QET (A,u)) in terms of the light-cone wave functions: 

F™(\,u) = ^f^dufidte^ [w P (u) - - ^(«))^k i0 *', (30) 

iT(A,fi>) = |^d M ^dte^[^H + (y^H-^2H)^ 



2 , ^ d 2 ,.A 



d-u 

e"^*, (31) 



where here and below u — 1 — u, and the superscript LO denotes the leading order in the 
1/rriQ expansion. 

Setting i^adr e q Ua i ^hqet dgf^gg the LC-HQET sum rules. The quark-hadron duality 
assumption is used to substitute the unknown resonances by the HQET result after a 
dispersion integration above some given threshold u> c , 



1 f°° Imb w ( A, v ) 
resonances = — du - — — + subtraction . (32) 

7T Juj c V — UJ 

The HQET spectral density can be obtained by the following double Borel transformation, 

1 ImF^ QET (A, u) = §W F« QET (A, Co) (33) 



7T 

with the Borel transformation being defined as 

B™ = lim Y t^L^l_ (34) 

Y x^oo, y= ^* flxed r(n) dX n v ; 

whose property, B^e pu} = 8(p — 1/t), is very useful in the calculation here. Evidently, 
before doing Borel transformations, one should first perform the Wick rotation on t. 

Finally, by performing a Borel transformation on both sides of the sum rule, so 

as to enhance the ground state contribution, suppress higher twist terms and remove the 
subtraction, one can obtain the sum rules, 

L a (X) = — e 2A / T / d v - ImF a LO (A, v) e^ T , (35) 
lb Jo IT 
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FIG. 2: The diagrammatic representation of the 2-point vacuum-pion correlation function Eq. ^dj), 
where the heavy solid lines represent the heavy quark, and the helical lines a gluon. 



where T is the Borel parameter, and 



r 2/v 



7T 2 
7T L 



^ImF a LO (A, v) = ^ 0(u Q ) ^(p p (u ) + v4(«o) ~ T^M 



A 



3A 



A 2 



3A 



A 2 



(36) 



(37) 



Here and in the following, m = 1 — v, and the prime denotes derivatives. After integration 
by parts, we obtain the final sum rules, 



L a {v-p) 



2F \ Jo 



/^ p («) j^-WW + -g 2 iu) 



-2uv-p/T 



- —92(0) 
o v-p 



-26vp/T 



(38) 



L b (vp) = ^z 2K/T {v-P J o 6 du 



<p v (u) + ^<fa(u) - ^gi(u) 2 — 

SI 1 z v ■ pi 



92{U) 



-2uv-p/T 



+ 



/ dgifj) 

6 r JM ' ' t>-p ^ dw 



92(d) 



-28v-p/T 



(39) 



where 6 = Min(l, and 6=1 — 6. We see that the LC sum rules become meaningless 
for processes with a very soft pion which would enhance higher twist contributions. 



B. 1/m.Q corrections from the effective current 

To obtain the 1/toq corrections from the effective current, we consider the 2-point 
vacuum-pion correlation function containing a covariant derivative, shown in Fig. ^ 

E?(\,u>) = 1 j c^xe^ (>K(p)\T{u(x)(T m iD)^h v (x), ^(0)i T5 d(0)} |0> 

= E^(X,Co)v lJ , + E^(\Co)p lx , (40) 
where m = 1, 3 is an index introduced for convenience, and 

(r 1 iD) M = 7 M 7^, (THD)^ = iD^. (41) 
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The hadronic representation of this 2-point correlation function can be expressed as 



£™ Hadr (A,£) = 2F " v 1 + resonances 
2A — uj 



(42) 



where the T™{X) have been defined in Eq. (pO|) . 

The above correlation function can be calculated directly in HQET. We choose the rele- 
vant momentum k to be parallel to v, i.e., = (Cj/2) v^, and work in the fixed-point gauge 
where A^(x) can be expressed as 



Ap(x) = x a I dw w Gnsiwx 



a/3 [ 



(43) 



Eq.(fiOD then can be given in terms of the above WJF(2)s, and the following 3-particle wave 
functions 0, [LjJ: 

(ir(p) \u(x)a fJil/ j 5 g s G a p(wx)d(0)\0) 



(7r(p) \u(x)j tl j 5 g s G a p(wx)d(0) |0) 



P» 



/ X a Pn\ 
\9afi J 



px 



Pa\9pii 



Xf3Py 

px 



VciiLp±(ai)e 



i{a\+waz) px 



+ U^ipaxp ~ Ppxa) / Von n (a i )e i ^ +wa ^ px . 
px J 

(ir(p) \u(x)j fl g s G a / 3 (wx)d(0)\0) 



if-* 



v vx ' v vx ' \ J 

i(oL\+waz) px 



+ ifn — {p a xp -Ppx a ) / Voii(p\\(ai)e l 
px J 



(44) 



where G a p = \e OL {i ar G ar , and Va^ = daida 2 da 3 5(l — a,\ — ct 2 — a 3 ). The wave function 
<^37r(tti) = (P3n(c<i, a 2 , a 3 ) is twist 3, and ip±, <f», y5j_ and <p» are twist 4. Using the identities, 
Eq. (tik and 



lulalp = 7/^a/3 ~ lag^p + ipdim ~ 1 £ m p5l 75 ■ 



(45) 



one obtains 



1 rOO 



^ du dte 1 ^ -H T7 



2 7o Jo 



d / Ar 



dt V 2 



^(m) + t 2 gi(u) p^ + fj^(p p (u)v^ 



A 



12' 



Hit 



+fi w ip p (u) + it[—\ip a (u) - g 2 (u) 



+ iX 



/ivr 

4 



^o-(w) - itg x {u) 



-2ig 2 (u)vu \ e 



V<\, ! dw ir j dte***-t 



- v^n (a,) k + 2^ ± ( Qti )^ } e- i(1 - Ql -^ 3) t* 



A ; y L V ? 37r(ai) 
Jit 



(46) 
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k 

2 Jo 



du / die*?* 



d / A r 



dt V 2 



_ A ^ 



A 



<p n (u) + t 2 g\{u) + ^ P {u) ) + \tg 1 (u)v i 



+ ig 2 {u){pu-v„) \e-^ -if /P a< / dww / dte^t 



A 



(47) 



where terms proportional to J? 2 are omitted in the chiral limit. 

After performing the corresponding Borel transformations, we obtain the sum rules for 
the effective current corrections: 



2A/T fu)e i 

= —— / dv-ImE™(X,v)e^ T 

It JO 7T 



(48) 



where 



-Im^(A^) 

7T 



7T 



71 



ye(w )[w (/^p( M o) + ^V CT ( M o) - t# 2 ( m o)) - y#2(«o) 



2 /""o 



A Jo 



dcti 



1-hr j 

6 

1-ai 1 

da 3 — 2 (p±(ax, 1 - ai - at 3 , a 3 ) 
a 3 



2 



6(w ) - uo( ^tt(mo) + /Wp( m o) - jg"(u ) ) + ^L(p a (u ) 



4 



2 /" n o 



+ T02(«o) - T^( M o) + x , 

A A A jo 



/a 



A^ 

l-ai 1 

cfa 3 — 



x(</?||(ai, 1 - «i - a 3 , a 3 ) - A^V 37r (ai, 1 - «i - a 3 ,a 3 ) ) 



2 

2 /•"<> 



A^ 



A^ 



4 // \ 
A A Jo 



da,\ 



l-ai 1 

da 3 -g y?N 1 - ai - a 3 , a 3 ) 



(49) 



(50) 



(51) 



-\mEl{X,u) 

7T 



^e(Mo) [ - M ( ^tt(Mo) + VnV P {uo) ~ T9l( U o) ) + Tfi^o) 



2 

4 2 f u ° 



1— ai 



da 3 — 2 V2||(ai, 1 - «i -a 3 , a 3 ) . (52) 



(15 



Note that one consequence of the equations of motion in Eq. ([H]), = —J 7 ^ * s explicitly 
satisfied by the sum rules. Combining the above sum rules with Eq. ( P2[ ) yields, 



F 6 (vp) 



— <p a (u) - — . 
6 U Jo 



da i 

x tp 3w (ax, 1 - ai - a 3 ,a 3 ) 



l-ai 1 

<ia3 — 

u-ai tt 3 



-2uv-p/T 



(53) 
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FIG. 3: The diagrammatic representation of the leading order vacuum-pion correlation functions 
with Okin and O mag insertions, (see Eq.(54)), where the heavy solid lines represent the heavy quark, 
and the helical lines a gluon. 



C. corrections from the effective Lagrangian 

Finally, we consider the 1/niQ corrections to the correlation function Eq. (|2!f ) from the 
sub- leading operators 0^ in and O mag in Eq. (f|). They are depicted in Fig. [5| One has 

F^\x, u>) = if d 4 xe lk - x (7r(p)\T{u{x)^h v {x), i f d*y &(y) , h v {0)i l5 d{0)\ |0) 



= F^\x,<j)v lt + F^\\,a)p lt , (54) 
where C± = 0^i n + C mag O mSug . By denoting 5L a /2rriQ as the 1/rriQ corrections to L a , namely 

SL a = K a + C ma g S a , (55) 
the hadronic representation of this correlator can be written as Jl] 



(^)Hadr IF 5F L a (X) 2F5L a (X) AF5AL a (X) 

F a (A, W ) = 2A _- + T " (2A _a )2 + resonances , (56) 

where 5F/2rriQ and 5A/2rriQ are the 1/rriQ corrections for F and A |Tj| [19]], respectively. 

We now calculate the sub-leading correlation function Eq. ( p^) in HQET, starting from 
the relation 

9% = WF(4)s + 0(a s )WJ r (2)s, (57) 

where the symbol = indicates that this equation holds at the level of matrix elements in 
terms of LC wave functions. This relation can be found by considering Eq.(E3) and the 
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equation of motion for the gluon fields, d a G a afi = —f abc A ba G c a p + g s h v vpT a h v . One can 
drop all the terms containing d@Ap or A^Ap, because contributions from wave functions 
with more than 3 particles and from order a s contributions are physically quite small |§. 
Therefore, we have 

F^\\Cj) = % J d 4 xe ikx i J d 4 y (n(p)\ t( u(x)^h v (*), h v (y) [ - d 2 + 2ig s A^ 

+ C nafS ^(r a pG af, ]h v (v), ^(0)z 75 rf(0) } |0) .(58) 

Here the higher-order terms in the heavy quark propagator may be included, as displayed 
in Fig. [5] (c). However, direct evaluation shows their contribution to be zero. Using the 
identities Eqs. ( PUD and (551) together with 

(59) 



i {9^a9/3u - 9^9au) + [?{iv9m ~ Ooa,9tf) ~ + i ^^59^ 'o >75 , (60) 

and after some algebraic manipulations, we obtain 



2 Jo Jo I dt 2 



A 



12 



X(f a (u) - g 2 {u) (pa - Va) 



A r 



_d 

dt 



2it 



12 



\<p a {u) - g 2 (u) (p„ - vj 



+ 2 At gi(u)p^ -iu\[ -[tip w (u) +t g x (u) Jp M + ^tif p {u)v^ 



— iuX 



12 



A^ ct (m) - g 2 (u) (p„ - u„) 



At£i(it) + 2^2 (tt) 



\i z g^u) - itg 2 {u) )p fi \e m ^ - f w I Va { I ds s I dww I dte 



A 



A 



x-tp fl <2 (f\\(ai) +<p±(ai) - i(l - a*! - wa 3 ) -t ip\\(a>i) \ e 



-i{\—a-±—swoiz)~t 



mag 



2<p ± (ati) - <p\\(ai) + 2(p ± (oti 



+ 2<f±(a i )v lJ 



(61) 



Having all necessary results at hand, we obtain the sum rules resulting from order 1/mq 
power corrections to the effective Lagrangian: 



6L a {\) 



7T 



+ 



1 

YF Jo 

where the sub-leading HQET spectral density functions are given by 



2 6A/T-5F La, (A) 



(62) 
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—ImFa" 1 ® ' (A, v) = ^Q{uo)\ u 2 ^(p' (uo) + ^(p"(uo)-^g2(uo 



6 



2u 



(63) 
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+ - / dai / da 3 — 



2 



6 



12 



^2(^0) - y^Oo) 



ln( 



a 3 



uq — a\ 



( l) <P\\(aci, 1 - «i - a 3 , a 3 ) 



(ai, 1 - ai - a 3 ,a 3 ) 



(64) 



It is interesting to note that the operator O mag in fact does not contribute to the HQET 
spectral density functions even at order 1/rriQ. 

IV. NUMERICAL ANALYSIS 



Let us now analyze the sum rules numerically. We use the pion wave functions collected 
and carefully discussed in Ref. [H, ITH]. They read, 



V? 37r (a;,/i) 



6uu jl + a 2 {n)-[5(u - uf - 1] + a A (^) — [21{u - uf - 14(u - u) 2 + 1] J 

1 + B 2 (n) l[3(u - u) 2 - 1] + B A (p)-[35(u - uf - 30{u - u) 2 + 3] , 

2 8 

Quull + C 2 (/m)-[B(u - u) 2 - 1] + C A (n)—[21(u - u) 4 - 14 (w - u) 2 + 1] 
12 8 

-S 2 (fi)u 2 u 2 + -e{fi)5 2 {fi)[uu{2 + 13m/) + 10m 3 lnu(2 - 3m + -u 2 ) 

2 2 5 

6 

+ 10w 3 lnw(2-3w + -w 2 )], 



—5 2 (fi)uu(u — w) , 
3 

360aia 2 a 2 1 + ^i,o(^)- (7a 3 — 3) + co> 2i0 (/i)(2 — 4«ia 2 — 8a 3 + 8a 2 ) 
+co , i i i(/i)(3aia2 — 2« 3 + 3a 2 



(p±{a h ii) = 30<5 2 O)(ai - a 2 )a 2 3 [- + 2e(/x)(l - 2a 3 )] 

(^||(aj,//) = 1205 2 (/i)£:(/i)(ai - a 2 )aia 2 a; 3 , 

^±(oi,/i) = 305 2 (/i)a 2 (l-a 3 )[^ + 2e(/i)(l-2a 3 )], 

<^ll(ai,/i) = -1205 2 aia 2 a 3 [^ + s(/j)(1 - 3a 3 )] , 



(65) 



where /x is the renormalization scale. The scale dependence of these wave functions is given 
by perturbative QCD. 

The light-cone sum rules obtained in HQET, and the pion wave functions in Eq. (|B"5D 
depend on the subtraction point fi. However, in the lowest order of a s , the value of the 
scale /i is ambiguous. Within HQET, a reasonable choice is ji ~ 2A. Thus we set the 
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renormalization point at /j,q = 1 GeV. In this (Ho) = 1.65 GeV, and r(/z ) 

The coefficients appearing in the above equation are then taken as H 



1.08. 



<h(lM)) = 0.44, a 4 (// ) = 
C^o) = 0.10, Gt(/i ) = 
/ 37 r(^o) = 0.0035 GeV 2 , c^ (/i ; 



= 0.25, ^(/io) = 0.48, 5 4 (/i ) = 1.15, 
0.067, 5 2 (/i ) = 0.2 GeV 2 , e(// ) = 0.5 , 
= -2.88, ^2,0(^0) = 10.5, = 0. 



(66) 



The pion decay constant / w = 0.132 GeV is taken from experiment. The parameters for 
heavy mesons are obtained from HQET sum rules [0| [L9J: A = 0.57 GeV, F(lGeV) = 
0.46 GeV 3/2 , SA = -0.35 GeV 2 and SF = -1.92 GeV. 



A. Universal functions 

In order to make the sum rules meaningful, we must first check the existence of sum rule 
windows, which are roughly given by Aqcd < T < 2A. The lower and upper limit will be 
obtained by the physical requirement that the Borel parameter T must be large enough to 
ensure that the higher-twist wave function contributions are suppressed, and at the same 
time small enough in order to make the resonance contribution not too large. To be specific, 
for the typical threshold values u c ~ 2 GeV, we find that, setting the Borel parameter to 
0.60 GeV < T < 1.00 GeV for the sum rules L a and L b , the twist-4 wave functions give 
contributions less than 17% and 3%, while the resonance contributions are lower than 27% 
and 6%, respectively. The two sum rules are quite stable in this region, as shown in Fig. |](a) 
and (b) where v ■ p is fixed at 2.0 GeV. Nevertheless, our numerical calculations show that 
the stability of the sum rules is not very sensitive to a change of u> c . In order to be consistent 
with HQET sum rule calculations, we adopt the same stability region of u c as that from the 
A sum rule |T^, [TJJ, i. e., u c = 2.0 ± 0.3 GeV. Fig. ^|(a) and (b) present the final results for 
L a (v ■ p) and L b {y ■ p), respectively, where the central value of the Borel parameter T = 0.80 
GeV is used. Note that the LC-HQET sum rules are meaningless in the soft pion region as 
mentioned above. To stay away from the soft pion region, we take v ■ p > 1.2 GeV. 

For the convenience of further applications, we parametrize the results by the following 
formulae in the region 1.2 GeV < v ■ p < 2.64 GeV, 

LT (v- p) = — l — - 2 , (67) 

a + a\v ■ p + a 2 (v ■ p) 2 

Lf\vp) = b + b lV -p + -^. (68) 

v ■ p 

Best fit values in which the maximal errors are less than 2% yields the following values of 
the parameters: 



UJ C [GeV] 


a [GeV- 1 / 2 ] 


ax [Gov- 3 / 2 ] 


a 2 [GeV- 5 / 2 ] 


b [GeV 1 - 72 ] 


bi [Gov- 1 / 2 ] 


b 2 [GeV 3 / 2 ] 


2.3 


-0.78 


4.00 


0.891 


0.403 


-0.0221 


-0.0385 


2.0 


-3.03 


6.71 


0.286 


0.409 


-0.0263 


-0.0448 


1.7 


-6.25 


10.9 


-0.743 


0.320 


-0.0049 


0.0343 



Next we consider the 1/ttiq corrections from the effective current. In the relevant universal 
function F 6 , only two twist-3 wave functions appear. Thus, in order to determine the 
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FIG. 4: Variation of the sum rule predictions on L a ^ with the Borel parameter T at various values 
of the continuum threshold parameter u c . 



sum rule window, especially its lower limit, one may consult the sum rules for T\ b that 
will contribute to the form factors directly while the renormalization-group effects are not 
considered. For threshold values of uj c ~ 1.8 GeV, requiring that both the twist-4 wave 
functions and the resonance contributions do not exceed 40%, the sum rule window can 
be determined to be 0.50 GeV < T < 0.70 GeV, in which we find that the resonance 
contributions for F% are less than 8%. The numerical results for the universal function Fq as 
a function of the Borel parameter T are shown in Fig. |6], where v ■ p is fixed to be 2.0 GeV. 
Fig. presents the results for F 6 (t> • p), where the central value T = 0.60 GeV is taken for 
the Borel parameter. 

The results can be parametrized by the formula below valid in the region 1.2 GeV < 
v ■ p < 2.64 GeV, with maximal errors less than 3%, 

F™(vp) = -. -. (69) 

6 v ; co + cif -p + c 2 (v-p) 2 v 1 
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FIG. 5: The universal functions L a ^(v ■ p) obtained from LC-HQET sum rules at various values 
of the continuum threshold parameter oj c . 

The parameters Cj are given by 



UJ C [GeV] 


C [GeV- 3 / 2 ] 


Ci [GcV- 5 /2] 


C 2 [GeV- 7 / 2 ] 


2.1 


-112 


30.5 


-4.93 


1.8 


-86.1 


1.54 


2.22 


1.5 


-57.6 


-34.0 


11.7 



We next analyze the sum rules resulting from the 1/mg insertions due to the power 
corrections to the effective Lagrangian. We find that, for the sum rule for 5L a , the twist-4 
wave function contributions are very small everywhere in the physically appropiate region 
of T, but the resonance contributions grow rapidly as T becomes large. Requiring the latter 
contributions to be less than 40% yields T < 0.55 GeV. Consequently, we set the range of 
the Borel parameter at 0.35 GeV < T < 0.55 GeV. The results of the sum rule for 5L b are 
very similar to the leading order sum rules. In the sum rule window 0.60 GeV <T< 1.00 
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FIG. 6: Dependence of the sum rule predictions of Fq on the Borel parameter T at various values 
of the continuum threshold parameter uj c . 
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FIG. 7: The sub-leading universal function Fq(v ■ p) obtained from the LC-HQET sum rules at 
various values of the continuum threshold parameter lo c . 

GeV and for threshold values u c ~ 1.8 GeV, the twist-4 wave function contributions are less 
than 11% and the resonance contribution does not exceed 18%. The stability of the sum 
rules for 5L a (a = a, b) with regard to variations of the Borel parameter T are shown in 
Fig. ||(a) and (b), where v -p is fixed at 2.0 GeV. Fig. ||(a) an d (b) present our final results 
for SL a (v -p) and SL^v -p), respectively, where the corresponding central value of the Borel 
parameter T = 0.45 and 0.80 GeV are used. It should be mentioned that, in these two 
sub-leading universal functions, the sub-leading HQET spectral density functions give small 
and negative contributions to the form factors, while the corrections coming from the SF 
and 5 A contributions are large, and positive when the sum is taken (for SL a , the proporation 
is about —1 : 2, and SL b about —1:5). 

We have parametrized these results by the following formulae in the region 1.2 GeV < 
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FIG. 8: Dependence of the sum rule predictions 5L a ^ on the Borel parameter T at various values 
of the continuum threshold parameter uj c . 



vp< 2.64 GeV, 

5Lf\vp) = 6 + M-P + 
The best fitting parameters are 



a + a x v ■ p + a 2 (v ■ p) 2 ' 
v ■ p 



(70) 
(71) 



U) c [GeV] 


a [GeV" 3 / 2 ] 


ai [Gcv- 5 / 2 ] 


a 2 [Gcv- 7 ' 2 ] 


b [GeV 3 / 2 ] 
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b 2 [GeV 5//2 ] 
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FIG. 9: The sub-leading universal functions 6L a &(u • p) obtained from LC-HQET sum rules at 
various values of the continuum threshold parameter lo c . 

where the maximal fitting errors are less than 5%. 

Finally we remark on the errors in the sum rules. The uncertainties of the continuum 
threshold parameter uj c induce errors (from both the LC-HQET sum rules themselves and 
the sum rules of A and F ||13|| ) less than ± 10% in the final sum rule results. However, one 
should keep in mind that the sum rule method typically has a ± (10 ~ 30)% uncertainty 
resulting from the duality assumption, uncertainties in the wave functions and other input 
parameters. 

B. Form factors to order I/tuq 

Having obtained the leading and the relevant sub-leading universal functions, we can 
construct the form factors fi } i(p ' p) from Eq. ([18]) for the B — > ir transition to order 1/nib- 
The results are shown in Fig. |10| (a) and (b), where we take m& = 4.7 GeV. For comparison, 
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the leading order results corresponding to u c = 2.0 GeV are also shown in this figure. Here 
the sub- leading universal functions roughly give a 11 ~ 23% and 5 ~ 12% enhancement for 
the form factors fi~* n {v ■ p) and /^"^(u • p), respectively. 
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FIG. 10: The form factors f^2* n (v ■ p) obtained from the LC-HQET sum rules at various values 
of the continuum threshold parameter uj c . 



V. COMPARISON WITH OTHER MODEL CALCULATIONS, DISCUSSION 
AND SUMMARY 



It is interesting to compare our results with those of other approaches. First let us look at 
the soft pion limit. The heavy meson chiral perturbation theory (HMxPT) [[20| (see also pT[ ) 
describes f\(v ■ p) and f^{v ■ p) in a single pole form in the soft pion region. Moreover, the 



20 



1/rriQ corrections do not change this behaviours 
they can be expressed by 



22]. In the case of B — > ir transitions, 



/i fl_Mr (v ■ p) 

J i HMvPT V J 



- HMxPT 



(v ■ p) 



CiF 

2/. 
CiF 

2U 



(1 + 



5F_ 

2m b 
SF 
2uih 



_ gv -p 

v ■ p + A B 

gvp 



P + A 



13 



(72) 
(73) 



where A# ~ mg. — mg = 0.045 GeV, g ~ 0.3 is the coupling of the pion to the heavy 
meson. It is easy to see that our result for fi~^ n (v-p) does not match to the HMxPT result. 
Concerning f^'^ip ■ p) our sum rule result Eq. (|39|) can be seen to reasonably well match 
on to the soft pion result considering the uncertainties in the values of F and g. By taking 
g = 0.28, one can find that the extrapolation of HM^P to large v -p matches quite well with 
the sum rule calculations at intermediate pion energies as shown in Fig [11] (a). For practical 
purposes, we use the following Gaussian-type function to make a smooth connection between 



the sum rule result of f± 



p) and that of HM^PT, 



fFc^v-p) 



C X F 



5F 



9o + flie 



2f„ V- ' 2m b 

( for 0.25 GeV < v ■ p < 1.2 GeV 



(74) 



where g = 0.125, g\ = 0.751, g 2 = 3.0 GeV -2 and 8 B = 0.23 GeV were obtained by 
matching both sides. This is plotted in Fig. |TT] (b). It should be mentioned that the 
matching results obtained here correspond to the case of central values of u c . For further 
support of our matching procedure we have included in the lattice NRQCD results from the 
JLQCD collaboration They can be seen to be well consistent with our matching result. 

The large pion energy limit is another interesting limit to consider because the large 
energy effective theory (LEET) [^4], ^] provides some model-independent information on 
this kinematical region. In the B — > 7r semileptonic decay, the light non-spectator quark gets 
a large amount of energy from the decay of the bottom quark. The light spectator quark 
system interacts with the energetic quark mainly at the energy scale Aqcd which is basically 
fixed by the size of hadrons. New symmetries appear in the limit of v -p/Aqcd ~~ * oo. They 



are subject to corrections by hard gluon exchange [ 26[ . In the leading order of the heavy 
quark expansion our results are 



oo 



0, L b (oo) = ^Tfa(l)I . 



where 



4 Jo 



(75) 



(76) 



This agrees with LEET. For the I/tuq corrections, our results are 
F fi foo) 



5L a (oo) = 0, 
5L h (oo) = 4r 



(77) 



T 



(25A/T-5F)I -Th) 0^(1). 
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FIG. 11: The form factors ff?*(v -p) as obtained from LC-HQET sum rules including a matching 
to the HMxPT result in the soft pion region. 



This can be compared with the LEET results |27J only after the hard gluon effects have been 
incorporated into these and can thus only been done in the future. There even is a school 
which assumes that the perturbative contribution is dominant for the B —>■ it transition [2S|. 
It is obvious that the LEET limit deserves further studies which certainly will be done in 
the near future. 

It is also interesting to compare our method and our results to the full QCD LC sum 
rule calculation. The authors of Ref. in fact considered the heavy quark limit of their 
full QCD LC sum rule calculation. Their results differ from L a (v ■ p) and L b {y ■ p) only by 
a simple transformation which were given in Eqs. (47) and (46) in their paper. There is a 
subtle difference from our results, though, which can be seen by letting %(1 — u)/u= v. It 
is only after taking vjv-p — > that the results of the two calculations fully agree. This limit 
may need further understanding. Numerically the QCD LC sum rules gave stable results for 
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q 2 < 17 GeV 2 which is consistent with the HQET LC sum rules requiring v ■ p > uj c /2 ~ 1 
GeV. These two methods could be in principle the same, provided that all the sub-leading 
corrections have been included in the calculations. Lacking such powerful calculations, the 
effective theory calculates a physical quantity in a most thorough and least complicated 
way through clearly separating the perturbative and nonperturbative parts of the quantity. 
Therefore in certain appropriate regions of the pion's phase space the HQET calculation 
may give more reliable results. Practically speaking the method is to calculate form factors 
at the hadronic scale, and the perturbative contribution is accounted for by multiplying in 
renormalization factors. 

It is physically useful to reconstruct the conventional form factors defined in Eq. ([[]) 
from the results combining LC-HQET sum rules and HJVllPT, by using the relation Eq.(§). 
In Fig. [12] (a) we present f+^iq 2 ), which is directly measurable in semi-leptonic decay 
involving light leptons. In Fig. [12] (b) we present our result for the scalar form factor 
fo^ n {q 2 )- The scalar form factor contributes to the decay B — > -nfv T and also enters the 
factorized amplitudes in non-leptonic two-body B decays. For comparison, the results from 



the lattice QCD simulations by APE J3tj, UKQCD f3T|, FNAL [g2j and JLQCD |23] are 
also shown in the figures. 

In Table | we compare our results to the results of other model calculations which includes 



the conventional HQET QCD sum rule calculation of Ref. [33]. We find that our results 



for f+^(q 2 ) are about 14 ~ 33% larger than those from the full QCD LC sum rule |], 



calculation which also includes a s corrections. 



TABLE I: Comparison between different approaches for the universal functions L a (v-p) and Lf>(v-p), 
and the form factor f+^ n (q 2 ) at some different kinematical points. 





L a (v-p) 


[GeV 1 / 2 ] 


L b {vp) [GeV 1 / 2 ] 








v ■ p [GeV] or q 2 [GeV 2 ] 


1.50 


2.64 


1.50 


2.64 


0.0 


6.0 


12.0 


LC-HQET SR (NLO) 










0.42 ± 0.04 


0.54 


0.77 


LC-HQET SR (LO) 


0.13 


0.06 


0.34 


0.32 


0.36 


0.47 


0.68 


HQET SR in LO || 


0.15 


0.13 


0.25 


0.21 


0.24 


0.32 




LC-QCD SR § H 










0.28 ±0.05 


0.42 


0.66 



Finally, in order to be concrete, we give our prediction for the decay width of the decay 



B c 



7T 



~e + v e using our form factors. We obtain 



V, 



ub 



4.08 x 10" 



1.46 ±0.30) x 10~ 16 GeV 



(78) 



By taking \V u b\ = (4.08 ± 1.18) x 10~ 3 from inclusive measurement of B — > X u £u |35 
we get T = (1.46 ± 0.30^72) x 10~ 16 GeV, where the uncertainties are from that of form 
factors and \V u b\ respectively. On the other hand, from the experimental result given in 



|, r B o = (1.55 ± 0.03) ps- 1 and Bi(B c 



7T e 1 v P 



+0.33 



24 ± 0.50) x 10 



-3 



1.8 ± 0.6) x 10~ 4 , we extract 



with the uncertainties from the form factors and the 



Ref. 

\V ub \ = (2.94 
experiments, respectively. 

Let us also comment on Ref. |ll| which gave leading order results for L a . In addition 
to that the analytical expressions of the sum rules differ from ours by a factor of 2, the 
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FIG. 12: The form factors /f "5 >7r (g 2 ) from combining LC-HQET sum rules and HM%PT, and from 
the lattice QCD simulations by APE |~~ 



UKQCD H, FNAL g|] and JLQCD |2J. 



numerical input used in Ref. [0 is quite different, such as the choice of the energy scale \i 
and the determination of the sum rule window. This naturally affects the final results on 
L a (vp). 

To summarize, we have applied the LC sum rule method to calculate the B iriu 
weak decay form factors to order I/tuq in the framework of HQET. We have calculated the 
leading and the relevant sub-leading universal form factors. Our form factor results have 
been matched to the appropiate soft pion results. We have also discussed the large pion 
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energy limit of our results. The full QCD LC sum rules and lattice QCD results have been 
compared with. In the future we are planning to include perturbative QCD corrections, 
try to incorporate the large energy effective theory into the QCD sum rule technique and 
perform detailed phenomenological analysis. 
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